Number of degrees of freedom of two-dimensional turbulence 

Chuong V. Trail and Luke Blackbourn 

School of Mathematics and Statistics, University of St Andrews, 
St Andrews KYI 6 9SS, United Kingdom 
(Dated: April 20, 2009) 

Abstract 

We derive upper bounds for the number of degrees of freedom of two-dimensional Navier-Stokes 
turbulence freely decaying from a smooth initial vorticity field uj{x,y,0) = cjq- This number, 
denoted by A^, is defined as the minimum dimension such that for n > N, arbitrary n-dimensional 
balls in phase space centred on the solution trajectory uj{x,y,t), for t > 0, contract under the 
dynamics of the system linearized about uj{x,y,t). In other words, N is the minimum number of 
greatest Lyapunov exponents whose sum becomes negative. It is found that N < CiRe when the 
phase space is endowed with the energy norm, and N < C2Rei^ + In iig)^/'^ when the phase space 
is endowed with the enstrophy norm. Here Ci and C2 are constant and Re is the Reynolds number 
defined in terms of ujq, the system length scale, and the viscosity z/. The linear (or nearly linear) 
dependence of N on R^ is consistent with the estimate for the number of active modes deduced from 
a recent mathematical bound for the viscous dissipation wave number. This result is in a sharp 
contrast to the forced case, for which well-known estimates for the Hausdorff dimension Dh of the 
global attractor scale highly superlinearly with We argue that the "extra" dependence of Dh 
on is not an intrinsic property of the turbulent dynamics. Rather, it is a "removable artifact," 
brought about by the use of a time-independent forcing as a model for energy and enstrophy 
injection that drives the turbulence. 

PACS numbers: 47.27.-1, 05.45.-a 
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I. INTRODUCTION 



Chaotic dynamics are characterized by the stretching and folding of volume elements 
in phase space (solution space). In the presence of dissipation, these can be accompanied 
by volume contraction. For a finite-dimensional system, volume elements can eventually 
collapse onto complex sets of zero volume having fractal structures, whose generalized di- 
mensions, such as the box-counting and Hausdorff dimensions, are significantly lower than 
the phase space dimension. For infinite-dimensional systems, volume contraction can occur 
for finite-dimensional volume elements. Furthermore, given a sufficiently large positive inte- 
ger (depending on physical parameters and initial conditions), this contraction can occur 
for arbitrarily oriented ra-dimensional volume elements following a trajectory — solution 
"curve" in function phase space — provided that n > N. This is the case if the sum of 
the largest N Lyapunov exponents at each point of the trajectory Ai > A2 > ■ ■ ■ > Aat, 
which can possibly be different for different trajectories, is negative. The smallest (which 
will be denoted by still) satisfying this condition thus defines the minimum dimension in 
phase space for which all n- dimensional {n > N) volume elements along a given trajectory 
contract during the course of evolution. This volume contraction means that the chaotic 
nature of the local dynamics can be "captured" and "contained" within a linear subspace 
having dimension not higher than A^. (This subspace may continuously change along the 
trajectory, though its dimension does not exceed A^.) For this reason, A^ can be thought 
of as an effective dimension of the dynamical system in question, in the sense that its local 
dynamics can be adequately described by an A^-dimensional model. When an attractor (or a 
global attractor) exists and A^ is common to every trajectory having initial data containing 
;he attractor, its box-counting and Hausdorff dimensions are both bounded from above by A^ 
l|, which is conveniently defined as the number of degrees of freedom. More precisely, these 
are bounded from above by the Lyapunov dimension D^, which satisfies A^ — 1 < < N 
and is defined by 

m 



In this study we determine upper bounds for A^ for two-dimensional Navier-Stokes turbu- 
lence freely decaying from a smooth initial vorticity field c^o in a doubly periodic domain of 
length scale L. Note that the global attractor for this case is trivial and has zero dimension. 
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However, the present problem is nontrivial because it is concerned with transient dynamics, 
most importantly during the stage of fully developed turbulence. The bounds obtained are 
expressible in terms of physical parameters and found to scale linearly or almost linearly 
(depending on the chosen norms for the phase space) with the Reynolds number R^, which 
is defined in terms of Uq, L, and the viscosity u. On the one hand, such scaling behaviors 
are in accord with heuristic arguments based on physical and mathematical estimates of 
the viscous dissipation wave number. On the other hand, these are in a sharp contrast to 
the forced case, for which well-known upper bounds for the Hausdorff dimension Dh of the 
global attractor have superlinear dependence on z/~ . We discuss this discrepancy 

and argue that the superlinear dependence of on z/~^ is not an intrinsic property of the 
turbulent dynamics. Rather, it appears to be a "removable artifact," brought about by the 
particular form of the forcing term used as a model for energy and enstrophy injection that 
drives the turbulence. Indeed, the "extra" dependence of on z/~^ would be removed if 
the energy and enstrophy injection could be made viscosity independent (more precisely if 
the injection could be bounded independently of viscosity), provided that this forcing model 
does not jeopardize the existence of the global attractor. 

II. PRELIMINARIES 

In this section we briefly recall a recently derived upper bound [7] for the enstrophy 
dissipation wave number k^. We deduce from this result an estimate for the number of 
active modes, by counting all modes having length scales larger than the dissipation length 
scale corresponding to kd- A brief functional setting of the two-dimensional Navier-Stokes 
system in the vorticity and stream function formulation is described, and the problem of 
phase space volume evolution is formulated. We avoid technical detail and use informal 
language. 

A. Number of active modes 

The two-dimensional Navier-Stokes system written in terms of the stream function ij) 
and vorticity uj = Aip is 

ujt + Ji'ip, uj) = uAu, (2) 
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where Ji^ip, uj) = ipx^^y — i'y^^x and u is the viscosity. We consider Eq. ([2]) in a doubly periodic 
domain of size 27tL. The initial vorticity field ujq is assumed to be smooth and have zero 
average. Equation ([2]) preserves the zero-mean property. This, together with periodicity, 
allows UJ (and tp) to be expressible as a Fourier series in terms of sin L~^{ix + my) and 
cos L~^{ix + my), where i and m are integers not simultaneously zero. In other words, 
the infinite-dimensional space (solution function space) can be spanned by the infinite basis 
{smL~^{ix + my), cos L~^{ix + my)}. Here {i,m) can be identified with a lattice of unit 
spacing on the upper half plane, with points on either half of the horizontal axis, including 
the origin, removed. 

The advection term J{il),uj) has many conservation laws. In particular, the total kinetic 
energy 11^11^2 = (|V^/'p)/2 = / dxdy /2, the total enstrophy ||u;f /2, and the peak 
vorticity are conserved. These are the most important conserved quantities and play 

prominent roles in the theory of turbulence. Under viscous effects, all these quantities 
decay, though in general at different rates. The enstrophy decays most rapidly, while the 
kinetic energy and the peak vorticity are far better conserved, with the latter probably best 
conserved j8|. 

For relatively small v, the free decay of a general smooth vorticity field presumably 
becomes turbulent, featuring a wide range of dynamically interacting scales that extend to 
the viscous dissipation range. This range is characterized by the dissipation wave number ky, 
which, according to the phenomenological theory of turbulence jd], is given by ky = x^^^/ v^^"^ ■ 
Here v = v ||Vc<j||^ /(47r^L^) denotes the enstrophy dissipation rate per unit area. Recently, 
Tran [71] derived the upper bound 

V ||Vu;||^ < , (3) 

for the dissipation rate i/UVcuH at its peak. Since both vorticity norms on the right-hand 
side of Eq. ([3]) decay we have the bound 

IIV-II < (4) 

which is valid uniformly in time, and the bound 

IIV7 II II l|l/2 

,, = ^<W^, (5) 

which is valid at least up to (and probably beyond) the time of peak enstrophy dissipation. 
The bound for the newly defined enstrophy dissipation wave number k^, compares favorably 



to ki, as it could be significantly smaller than k^, By the very definition ([5j), enstrophy 
dissipation is strongest in the vicinity of kd- The wave numbers greater than k^ are effectively 
suppressed by viscous forces and virtually inactive. The number of dynamically active modes 
Nc corresponding to k < k^, are therefore given by 



Nc^ji< (6) 



where ko = 1/L is the smallest wave number. The quantity L Hcuolloo ^^^Y be identified with 
the fiuid velocity. Perhaps, \\uJo\\ is a better representative of the fiuid velocity; nevertheless, 
when it comes to the definition of the Reynolds number R^, we use L \\ujq\\^ and \\ujq\\ 
interchangeably. With this identification, the term on the right-hand side of Eq. (jS]) may 
be defined as the Reynolds number Rg. Hence, Eq. ([6]) can be rewritten in a more compact 
form 

iVc < Re. (7) 

From our experience in numerical simulations of two-dimensional turbulence, the estimate 
([7]) is sharp — in fact, spot on. For example, for the standard numerical domain 2tt x 27r 
and an initial vorticity maximum Hcuolloo ~ ^tt, the simulations of Dritschel, Tran, and Scott 
8] using 47r(8/3)^/z/ grid points adequately resolve the dissipation scales. This resolution is 
obviously consistent with Eq. ([71), within an order of magnitude. As will be seen in the next 
section, the estimate (171) for Nc fully agrees with the number of degrees of freedom discussed 
above. 



B. Problem formulation 

The problem of phase space volume element contraction (or expansion) is intimately 
related to the stability of solution with respect to disturbances. To investigate this problem, 
we consider the linear evolution of a deviation (p of the stream function ip (corresponding to 
a deviation A0 of the vorticity u) governed by the linearised equation 

A<Pt + J{<P, uj) + J{iIj, A0) = V, (8) 

where uj (and ip) solves Eq. ([2]) with initial vorticity ujq (and initial stream function ipo). By 
taking the scalar product ((■)) of Eq. ([8]) with and A0 we obtain the respective evolution 
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equations for the energy norm ||V0|| and enstrophy norm ||A0||, 

l|V0|||||V0|| = (MV',A0))-i.||A0f (9) 

and 

II A0II I IIA0II = -(AM'/', ^)) - y II VA(/.f . (10) 

The respective exponential growth (or decay) rates A and A for ||V0|| and ||A0|| can be 
readily deduced and are given by 

A = I In II V0II = ((0J(V', m - V II A</,f ) (11) 

at IIV0II 

and 

A = I In II A0II = ^ ((AM0, ^)) + ^ II VA0f ) . (12) 

These rates provide a comprehensive picture of solution stability, quantitatively describing 
how solutions with nearby initial data disperse from one another. 

Two natural norms for the present problem are the energy and enstrophy norms. We 
will refer to the phase space equipped with the energy (enstrophy) norm as the energy 
(enstrophy) space. In the course of evolution, consider a trajectory commencing from a 
given initial condition. At an arbitrary point on the trajectory (i.e., at an arbitrary instance 
in time i > 0), we calculate the greatest growth rate A (A) and identify the corresponding 
most unstable "direction" by considering the problem of maximizing A (A) with respect to 
all admissible 0. We denote by (Ai,v9i) [(Ai,-;?!)] the solution of this problem, where for 
convenience {&) has been normahzed, i.e., ||V(^i|| = 1 (lA^j^iH = 1). The second greatest 
rate A2 (A2) and the corresponding second most unstable direction (/?2 {^Qt) orthogonal to </7i 
(■j?i) is obtained by the same maximization problem subject to the orthogonality constraint, 
i.e., (V(/?i- V(/?2) — ((Ai?iAi?2) ==0). By repeating this procedure n times, we obtain the set 
(/92, ■ ■ ■ , V^n} ({''^i) ■ ■ ■ ; "^^n}) of mutually oithonoimal functions and the corresponding 
set of ordered rates Ai > A2 > • ■ ■ > A„ (Ai > A2 > • • • > A„). These may be defined as the 
first n local Lyapunov exponents, and their existence is guaranteed since the maximization 
problems are expected to return unique solutions. Note that for the conventional Lyapunov 
exponents, existence can be a major issue, even for low-dimensional systems of a few degrees 
of freedom. 
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Now, in the linear subspace spanned by {ipi, ip2, ■ ■ ■ , fn} {{'^i, '^2, - ' ' ^ ^n}), consider an n- 
dimensional ball B{-,r) of radius r centred at the point discussed above. The n-dimensional 
volumes v (in the energy subspace) and V (in the enstrophy subspace) oi B{-,r) are given by 
f oc r" llVv^ill ||Vv22|| ■ ■ ■ llVc^rill = r'^ and \^ oc r" HA-i^iH ||A7?2|| ■ ■ • IIAi^nH = r", respectively. 
(See the book of Temam 10|] for a formal definition of volume based on the related concept 
of exterior product.) The respective equations governing the evolution of v and V under the 
linearised dynamics described by Eq. ([8]) are 

, n n 

J^^rlv = J2^^ = Y. ((^^'/(^, A^O) - ^ W^vf) (13) 

i=l i=l 

and 

, n n 

-ln\/ = ^A, = (^{A^'J{^\ u)) + V \VM'f^ . (14) 

j=l i=\ 



In deriving Eqs. f[T3|) and flMj) . we have used Eqs. f[TT|) and f[T2|) . respectively. The sum 
Sr=i -^i (SILi represents the exponential growth or decay rate of v (\^). When this sum 
is negative, the volume of the n-dimensional ball B{-,r) contracts exponentially. Note that 
by construction, B{-,r) is optimally "oriented" to be least contracting. This means that 
if Sr=i (Sr=i ^i) i^ negative, then volume contraction becomes universal for all n- or 
higher-dimensional balls locally centred at the point in question. Furthermore, if this point 
is taken arbitrarily on the trajectory, which is the case in this study, then volume contraction 
becomes universal along the trajectory. 

The determination of then reduces to minimizing n such that the sum on the right- 
hand side of Eqs. ( |T3l) and (HM is negative. We use the mathematical techniques developed 
in the 1980s by Babin and Vishik [4] and Constantin, Foias, and Temam [s], ^ 11] for 
estimating the attractor dimension of forced two-dimensional Navier-Stokes turbulence. See 
also .he pape. of Doen.g a.d G.bbon Q for the sa.e ..eat.ent the st.ea. function 
and vorticity setting. As can be seen in the next section, the derivation of upper bounds for 
A^ is equivalent to the determination of the Hausdorff dimension of the global attractor in 
the forced case. The main difference is that although the present formulation is specifically 
designed to handle the decaying case, which has a trivial global attractor, its scope of 
application is broad. In general, the present notion of degrees of freedom makes sense for 
general dissipative dynamical systems, provided that bounded solutions exist. There are 
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virtually no other technical requirements for the application of the method. In particular, 
no a priori knowledge of the existence of an attractor is required. 



III. RESULTS 



This section presents the calculations described above, leading to upper bounds for A^. 
The treatment is relatively self-contained. However, the reader, who is interested in further 
detail related to the analytic inequalities employed in various stages of the calculations, is 
referred to the cited papers and references therein. 



A. Degrees of freedom in energy space 



i=l 



We begin by deriving an upper bound for N in the energy space. From Eq. (fT3|) we have 

n n 
■ ' i=l 

n 

i=\ 
n 

i=\ 

n 

2=1 

n 

< E ((I VKI^xllV^.I + |^;||V^,|)) - u \\Ay,f 

i=l 

n 



i=l 



< M 



i=l 



i=l 



(15) 



where integration by parts and the Cauchy-Schwarz inequality have been used. For further 
estimates of the terms on the right-hand side of Eq. f|T5l) . we employ the following two 
analytic inequalities concerning the orthonormal set |g9^|"_i with respect to the energy 
norm. First, we have the Lieb-Thirring inequality 0, llo| 

1/2 



i=l 



(16) 



j=i 



where ci is a non-dimensional constant independent of the set {y?*}"^^. Second, we know that 
for n ^ 1, there are approximately n basis functions (trigonometric functions mentioned 
earlier) within the wave number radius y/n/L. Their (repeated) eigenvalues under —A 
are (£^ + m^)/L^, where i"^ + < n. These constitute the first n eigenvalues (in non- 
decreasing order) of —A and sum up to approximately n'^/L'^. It follows from the Rayleigh- 
Ritz principle that 

n 2 
i=l 

where C2 is another non-dimensional constant independent of the set {ip''}^^i. By substituting 
Eqs. ([16]) and into Eq. we obtain 

n / n \ 1/2 

Y.^^<\J2\\^V>f] (c^M-u'j-n). (18) 

i=l \ i=l J 

It follows that Yl^=i -^j ^ when n > ciL \\u!\\ /(c2Z/). Hence we deduce the bound 

N < < Ci^^ = C,R^, (19) 

where C\ = Ci/c2 and Re has been redefined by replacing L \\ujq\\^ with \\uJo\\. Note that the 
precise result should be that is no greater than the least integral upper bound for CiRe] 
however, in writing Eq. ( fT9l) . we have opted to ignore this exceedingly minor detail. Equation 
f|T9|) gives a clear linear dependence of N on R^.. Thus, we have essentially recovered the 
bound ([7]), up to the constant factor Ci and a slight difference in the definition of Re, which 
was obtained earlier by counting the active modes from the smallest wave number ko = 1/L 
to the dissipation wave number = HVcuH / ||^^||• 



B. Degrees of freedom in enstrophy space 

An upper bound for in the enstrophy space is derived in a similar manner. From Eq. 
iHM we have 

n n 

^A, = - ^[{M'j{&,uj)) + V ||VA^9'||') 

i=l i=l 

n 

< (^{\M'\\V^'\\Voj\) - z/||VA#^||^) 

i=l 
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1/2 



,i=l 



i=l 



i=l 



2\ 1/4 



(|Vu;|4/3)3/4_^,^||VA^?^| 



,i=l 



i=l 



< 



1/2 



i=l 



1/2 



Vuj\\ — V 



(20) 



where Holder's inequalities with the pairs of conjugate exponents 4/3 and 4 and 3/2 and 
3 have been used in the penultimate and final steps, respectively. For further estimates 
of the terms in this equation, we employ a few more analytic inequalities concerning the 
orthonormal set in the enstrophy space. First, we have 0, [ll| 



i=\ 



<cl l + lnE^^||VA^9' 



(21) 



i=l 



where C3 is a non-dimensional constant independent of the set {"(9*}"^^. Second, a version of 
Eq. f|T6|) for the present orthonormal set is 

1/2 



2 = 1 



<ci Ell^^^* 



(22) 



, j=i 



Finally, a version of Eq. f|T7|) for {i?*}"^^ is 



n 2 

^||VA?9'||'> 

i=l 



L2 



Now by substituting Eqs. (gD, (|2T]), and ([22]) into Eq. ([20]) we obtain 



(23) 



1/2 



C'(l + lnO^/ 



1/4 



l|l/2 II I 

^oIL ll^ol 



,i=l 



1/2 



'EllVA^r 



2=1 



^^^^^ /"r./ /1 , Ikolli^^ Ikol 



L2 



3/2 



(24) 



where C = v^CiCg, ^ = ELi l|VA^?^f , and = {L^ \\ujo\\^ \\uJo\fy/yu. Note that by 
Eq. fl23l) we have ^ > c^n^. Hence without the logarithmic term, it would be straightforward 
to substitute this into Eq. and deduce an upper bound for similar to Eq. (fT9!) with 
the newly defined Re replacing its previously defined (and comparable) counterpart. Since 
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we are interested in the case ^ ^ 1, the logarithmic term should introduce a small departure 
to the linear dependence of N on Re only. In order to account for In^, we can "cover" it by 
a fraction of ^, say ^/2. By elementary calculus, we find that 

C'Rl/^ (1 + InO'/' - ^ < V2C'Rl/^ (1 + InRef^ , (25) 



where we have dropped a negative term on the right-hand side. It follows that 
C'Rf^ (1 + InO'/' - < V2C'Rl/^ (1 + InRef^ - ^ 



< V2C'R'J' (1 + \nRef' - (26) 

The condition XliLi ^ is satisfied when the right-hand side of Eq. fl2Bl) is non-positive. 
This requires a straightforward condition for n which in turn yields the result 

N <C2Re{l + lnRe)^^\ (27) 

where C2 = {8C'^y/yc2. 

As expected, Eq. (!27|) gives an essentially linear scaling of with Re since the superlinear 
dependence on Re, due to the logarithmic term, is slight for large Re- Given that the same 
linear scaling was found earlier in the energy space, this is somewhat surprising. The reason is 
that the energy in two-dimensional turbulence is predominantly transferred to smaller wave 
numbers while the enstrophy is predominantly transferred to larger wave numbers. This 
undoubtedly implies that the enstrophy dynamics have relatively more degrees of freedom 
than the energy dynamics. Hence, it is somewhat counter-intuitive that Eqs. f|T9|) and fl27|) 
do not differ by much. A possible explanation is that Eq. (fT9|) may not be as optimal as Eq. 
f l27|l . Some qualitative support for this possibility turns up in the next subsection. 

C. Discussion 

In the 1980s, estimates were derived for the Hausdorff dimension Dh of the global at- 
Tactor of the two-dimensional Navier-Stokes system driven by a time-independent force / 
. These estimates have been known to be sharp, allowing just minor improvements 
for the attractor dimension in the energy space only [l^, ll^- present notations, the 

respective bounds for Dh in the energy and enstrophy spaces are given by 

11 



and 



C 



"6-2/3(1 + In G')^/^ (29) 



where c' and c" are constant and G is known as the generahsed Grasshof number. Although 
G has some certain physical significance, its highly superlinear dependence on z/~^ appears 
to make Eqs. (1281) and (1291) in disagreement with the bounds for Nc and for derived earlier. 
We claim that this apparent disagreement is due entirely to the particular form of / and 



Dorate on this claim. 
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could be fully reconciled. For the remainder of this paper, we wil 

Due to rigor requirement in the mathematical formulation 
forcing / has been used as a model for energy and enstrophy injection. The forced Navier- 
Stokes equations 



ut + {u- V)u + Vp = v^u + / (30) 
V • u = 



then admit the following evolution equations 



-'- II l|2 Win l|2 , / j?\ 

--HI =-H|Vn|| +{u.f) 

< -z/||Vwf + ||Vn|| ||V"^/|| 

<-^l|V.f + ^ (31) 



and 

1 d 



< -u\\Au\f + \\Au\\ ll/l 



|2 

2 ■ 2z/" 

|2 /n 1 ^1 iivv Il2 



< IIAnf + % (32) 



for the energy \\u\\ /2 and enstrophy ||Vw|| /2, respectively. In Eqs. fl3Tl) and fl32l) . the 

1 2 2 

terms ||V~ /|| /(2z/) and ||/|| /(2z/) represent upper bounds for the energy and enstrophy 
injection rates, respectively. Their dependence on u is inescapable because the injection 
rates {u ■ f) and —{Au ■ f) themselves are fiow dependent. 

We now demonstrate how the viscosity dependence of the injection rates (or more pre- 
cisely of the upper bounds for the injection rates) contributes to the superlinear scaling of 
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Dh with z/ ^. To this end, let us recall the intermediate steps , Q toward ([28]) and 
fl29|) given below 

^- -2V2 

I^^/ < c ^ ^ 33 

and 

rl/2\ 2/3 / ^1/2 \ 



^^<^ I " ^" I h + " ^" I ' (34) 

where the overline denotes the supremum of an asymptotic average. From Eqs. (l3T|l and 
we can deduce the forced dissipative balance equations 



< ^ (35) 



and 



■1/2 ^ 11/ 



||Au|r < (36) 

Upon substituting these into Eqs. (!33|) and we recover Eqs. (!28|) and (!29|) . respectively. 

However, if the driving force could somehow be modelled in such a way that the averaged 
energy and enstrophy injection rates would be bounded independently of viscosity, say by 

and rf, respectively, then Eqs. fl35l) and fl36l) would become 

IIV^II < ^ (37) 

and 



Ai^f <-^. (38) 



1/2 r] 



V 



1/2- 



Then upon substituting these into Eqs. ( I33l) and ( l34l) . we obtain 



< ^% (39) 



and 



(J- 2 \ 2/3 / r2 \ 1/3 

One can see that Eq. (HO!) has the desired scaling, i.e., linear dependence on v^^ with a 
logarithmic "correction" as in Eq. (p7|) . Hence, for the dimension estimate in the enstrophy 
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space, the "extra" dependence on z/~^ would be completely removed if the enstrophy injection 
could be bounded independently of viscosity. Note, however, that only part of the extra 
dependence on u^^ would be removed from the dimension estimate in the energy space 
when the energy injection is made independent of viscosity. This strengthens our earlier 
suggestion that the estimate for (and Dh) in the energy space might not be as optimal 
as its counterpart in the enstrophy space. 



Closely related to the present notion of degrees o: 
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freedom are the concepts of determining 
[isl . Q. A large body of research on 
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modes, nodes, and finite-volume elements 
the number of degrees of freedom deduced from these concepts has produced upper bounds 
proportional to G (cf. |l9|). Like the estimate for Dh in the energy space, these bounds would 
reduce to oc z/~'^/^ if the enstrophy injection could be bounded independently of viscosity. 



In another mathematical study 



22l | relevant to the present problem, it has been found that 



when the time-independent / in (l30l) consists of a single Fourier mode, the unstable manifold 
emanating from the stationary solution —A^^f/u has a dimension not lower than oc G"^^^. 
This is a lower bound for the Hausdorff dimension Dh of the global attractor, and the 
superlinear dependence on of this bound is a consequence of the dependence on u of the 
stationary solution. 

In passing, it is worth mentioning that numerical simulations of two-dimensional turbu- 
lence have routinely used a variety of forcing that provides steady energy and enstrophy 
injection rates and r/^. This class of forcing includes white noise and flow dependent 



forcing 
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211 ]. While such a class of forcing is numerically desirable and realistic in some 



sense, it may render Eq. ( ]30i) incompatible with the mathematical formulation leading to 
the desired estimate (HOj) . Nevertheless, for the present approach, there are no technical dif- 
ficulties in arriving at this estimate as an upper bound for the number of degrees of freedom 
in the present sense. 



IV. CONCLUSION 



In conclusion, we have derived upper bounds for the number of degrees of freedom N 
of two-dimensional Navier-Stokes turbulence freely evolving from a smooth initial vorticity 
field in a doubly periodic domain. This number is defined as the minimum dimension such 
that arbitrary phase space volume elements of no lower dimensions along the solution curve 
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in phase space contract exponentially under the linearized dynamics. This means that 
the (locally in time) turbulent dynamics could be sufficiently "contained" within a linear 
subspace whose dimension does not exceed N. In essence, N represents a reduced dimension 
that a modelled system should achieve in order to describe the turbulence adequately. It 
is found that < CiRg in the energy space and N < C2-Re(l + lni?e)^^^ in the enstrophy 
space. Here Ci and C2 are constant and Re is the Reynolds number, which is defined in 
terms of the initial vorticity, the system size, and the viscosity. These results are consistent 
with the number of active modes deduced from a recent mathematical estimate of the viscous 
dissipation wave number kd = ||Vc<j|| / ||c<j||. 

The present estimates for N have been compared with well-known bounds for the Haus- 
dorff dimension of the global attractor in the forced case, and the apparent difi^erence 
between the linear (or nearly so) scaling of N with Re and the highly superlinear dependence 
of Dh on the inverse viscosity has been discussed. We have argued that the superhnear 
dependence of Dh on i/'^ is not an intrinsic property of the turbulent dynamics and further 
suggested that this is a "removable artifact," arising from the use of a time-independent 
forcing as a model for energy and enstrophy injection that drives the turbulence. This sug- 
gestion has been strengthened by the fact that the "extra" dependence of on would 
be completely removed (at least for the estimate of Dh in the enstrophy space) if one could 
model the driving force in such a way that the enstrophy injection rate does not depend 
on the viscosity. Such a forcing can be seen to be more realistic than ones with viscosity 
dependent input. 

In the present analysis, we simply follow a trajectory starting from an arbitrary smooth 
initial vorticity field in the solution (function) space of the two-dimensional Navier-Stokes 
equations and monitor the evolution (under the linearized dynamics) of the volumes of n- 
dimensional balls centred on the trajectory. We estimate how large n should be to ensure 
that these volumes contract exponentially. This turns out to be equivalent to the method 
of estimating the Hausdorff dimension of the global attractor of the forced system. The 
present approach can be seen to be highly flexible in application. In general, it is applicable to 
either autonomous or non-autonomous, forced or unforced, and finite-dimensional or infinite- 
dimensional systems. There are virtually no special requirements, other than existence of 
solution, for the present definition (and method of analysis) of the number of degrees of 
freedom to make sense. In particular, the existence of the usual Lyapunov exponents is not 
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an issue. Furthermore, there is no need for a priori knowledge of the existence of an attractor 
(or a global attractor), whose generalized dimensions would normally be considered as the 
number of degrees of freedom of the dynamical system in question. Given all this, we may 
apply the present approach to less idealized and more realistic dynamical models without 
risking to compromise mathematical rigor. 

We thank the staff of the Isaac Newton Institute for Mathematical Sciences for their 
hospitality and support during the program on "The Nuture of High Reynolds Number 
Turbulence," when this paper was completed. 
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